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It is shown that the noise process in quantum computation can be described by spatially corre- 
lated decoherence and dissipation. We demonstrate that the conventional quantum error correcting 

QO ■ codes correcting for single-qubit errors are applicable for reducing spatially correlated noise. 
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Quantum computers hold the promise for solving many hard problems much more effectively than their classical 
counterparts [1-2]. However, in practice the inevitable noise and decoherence process will diminish the advantages of 
quantum computation [4]. To overcome the fragility of quantum information, many kinds of techniques have been 
discovered to combat noise and decoherence in quantum computers [5-18]. Among these techniques, quantum error 
correction is the most important one [5-12]. 

In quantum error correction schemes, the input state is encoded into a state in a larger Hilbcrt space so that 
it can be recovered from a certain class of errors. The error operators are identified from physical consideration. 
Conventionally, it is assumed that qubits (quantum bits) are decohered independently In this circumstance, after a 
short time interval, the most important errors are those described by single-qubit operators [19]. Hence, quantum 
error correcting codes (QECCs) need only to correct single-qubit or few-qubit errors [5-12]. 

The assumption of independent decoherence plays an important role in quantum error correction schemes [5,10]. 
Except independent decoherence, there are possibly other kinds of decoherence, for example, collective decoherence 
[15,16]. Alternate schemes, called quantum error avoiding codes, have been proposed for reducing collective decoher- 
ence [16-18]. Are quantum error correction techniques applicable for reducing collective decoherence? The answer is 
shown to be positive for a special collective dephasing model [20]. In this paper, we consider the most general type 
of noise of many qubits. It is shown that to a good approximation, the general noise process can be described by 
spatially correlated decoherence and dissipation, which includes independent decoherence and collective decoherence 
as its special cases. The master equation is derived, and we identify the error operators for general spatially corre- 
lated noise, using the quantum trajectory approach [21,22]. The error operators for spatially correlated noise arc no 
longer single-qubit operators. They are either a sum of single-qubit operators (the quantum jump errors) or a tensor 
product of them ( the effective evolution error). However, we demonstrate that the conventional QECCs correcting 
for single-qubit errors remain valid for reducing this kind of spatially correlated noise. 

We consider noise and decoherence of L qubits. The I qubit is described by the Pauli operator The most 
general noise process of L qubits can be described by the following interaction Hamiltonian (setting % = I): 

Hj (t) = 91 £ <r? (t) + .92 E E < ! ® °7 T 7 1 ' (*) 
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where ai — 1,2,3, and I = 1,2, ■■■,L (of' 2 ' 3 correspond to <rf' y ' z , respectively). All the Tf (t), T^T 1 ' (t) , and 
^i2° 2 N aN (*)> generally dependent of time, are noise terms, which may be classical stochastic variables or stochastic 
quantum operators, corresponding to classical noise or quantum noise, respectively. The Hamiltonian (1) includes all 
possible interaction terms between the qubits and the noisy environment, and it is very complicate. Fortunately, in 
practice this complicate description is not necessary. The coupling coefficients 51, 32, ■ ■ ■ and g n normally satisfy the 
condition g\ 3> 32 3> • • • 3> Ql- Hence, the most important noise always comes from the first term of the right hand 
side of Eq. (1), if this term does not reduce to zero due to some special symmetry. As a good approximation, we can 
safely drop off all the high order nonlinearities in Eq. (1). The interaction Hamiltonian is then simplified to 

ffj(t) = ff i5>?n"(t). (2) 

l, a 

The explicit expressions for the noise terms Tf (i) depend on the concrete physical model of quantum computation. 
However, it is reasonable to assume that Tf (t) satisfy the following conditions: 

<r? «)>,„, = 0, (3) 

where (} env denotes average over the environment. In Eq. (4), all the f"^ (t — t^j are correlation functions. From 
the Hcrmiticity of Tf (t) , we have 



=[/A Q (-*)J • ( 5 ) 
The general form of the master equation with the interaction Hamiltonian Hi (t) is expressed as [23] 

d 

—p (t) = - J drtr env {[Hj (t) , [Hj (t-r),p (t) ® p env ]]} , (6) 

where p env is the density operator of the environment, and p(t) denotes the reduced density operator of 

the qubits in the interaction picture. Substituting the Hamiltonian (2) into Eq. (6), and using the 

conditions (3) and (4), we derive the following master equaiton for noise and decoherencc of L qubits: 
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where the coefficients 
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From Eq. (5), it follows that the matrixes A = 



A i'i 



arc Hcrmitian. The first term of the right hand 



and B = 

side of Eq. (7) represents the environment-induced Lamb phase shift, and the second term represents decoherence 
and dissipation of the qubits. If the coefficients and B^ arc directly probational to S u > , Eq. (7) describes 
independent decoherence (in the terminology of Refs. [15-19]). In contrast, if the correlation terms A 1 ^" and Bf,° with 
I 7^ I attain the maximum, Eq. (7) represents collective decoherence [15,16]. In general circumstances, the qubits are 
subject to spatially correlated decoherence. 

We are interested in the problem that to what extent the noise described by Eq. (7) can be reduced by the 
conventional quantum error correction techniques. To examine the problem, we need first to identify all the first- 
order error operators for spatially correlated noise. It is convenient to use the quantum trajectory approach to 
attain this goal. The quantum trajectory approach is a recently-developed numerical simulation method for solving 
complicate open quantum systems [21,22]. In this approach, the evolution of the dissipative system is represented by 
an ensemble of wave functions that propagate according to the effective Hamiltonian interrupted by random quantum 
jumps [24]. To use the language of quantum trajectories, we need to re-express the master equation (7) in a diagonal 
form. The correlation matrix A is Hcrmitian, hence it can diagonalized by a unitary matrix U — [U^"], i.e., we have 



A l'l -2s U l'k ^k U kl i 

fc,7 



(10) 



where with k = 1,2, • • • ,L and 7 = 1,2,3, are eigenvalues of the positive-definite Hermitian matrix A, which 
should be positive real numbers. Define the operators sj. by the equation 



-2 =E W ■ 



(ii) 



In general, are no longer Hermitian operators. With the transformations (10) and (11), the master equation (7) is 
rewritten as 

|p (t) = -iH eff p (t) + (t) Hl ff + J2 &lp it) (stf , (12) 

fc,7 

where the non-Hermitian effective Hamiltonian is 

l,l' ,a,0 k '1 

The first term of the effective Hamiltonian is the Hermitian Lamb phase shift, which in general cannot be simplified 
by introducing the operators s^, since the matrixes A and B do not necessarily commute with each other. The second 
term of the right hand side of Eq. (13) is the non- Hermitian damping Hamiltonian. 

Suppose that in a finite time To, we perform N times error corrections. In a short time interval At = j^-, we need to 
find all the error operators up to the first order of At. In the language of quantum trajectories, the system evolution 
described by Eq. (12) is represented by an ensemble of pure states that evolve according to the effective Hamiltonian 
(13), interrupted at random times by quantum jumps sj.. Up to the first order of At, the normalized state after At 
will be cither 



|*Z(At)> = ^^*Z|*(0)) (14) 

with probability pi = ^ At (\I> (0)| (si) |<J> (0)) in case of a jump in decay channel (k, 7) at a random time in the 
interval At, or 

|*o (At)) = -L exp (-iH eff At) |* (0)) 



fc.7 1,1', a,0 



(15) 
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with probability Po = 1- E Zl At (°)l ( S D s l I* (°)> _ o (At 2 ) if no jump occurred. Let Q 
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?3(fc-i)+ 7 = \/ l^ s fc> and P3(fc-i)+ 7 = Pfc, where fc = 1, 2, • • • , £ 



and 7=1,2, 3. With this notation, the system state after a short time At is then represented by the following density 
operator 



3L 



p(At)^J2p^ lP (0)Qi + o(At 2 ), (16) 



n=0 



where p(0) = |\& (0)) (0)|. In the above equation, Q represents the effective evolution error, and Qk 
(k = 1, 2, • • • , 3L) represent the quantum jump errors. All the Q n (n = 0, 1, ■ ■ ■ , 3L) make a complete set of the 
hrst-ordcr error operators. 

For independent decoherencc, the correlation coefficients A^" and B^" are directly proportional to S u > . All the 
first-order errors Q n then reduce to single-qubit operators. For general spatially correlated decoherence, the first-order 
errors are no longer single-qubit operators. The quantum jump errors Qk (k = 1, 2, ■ ■ ■ , 3L) are expressed as sums of 
single-qubit operators, and more seriously, the effective evolution error Q includes the terms that are tensor products 
of single-qubit operators. Hence, it is not clear that this kind of decoherence can be reduced by the conventional 
QECCs. In fact, it has been suggested that to combat the effective evolution error, more involved and less efficient 
QECCs need be devised [25]. However, here we show that the conventional QECCs correcting for single-qubit errors 
remain applicable for reducing general spatially correlated decoherence, if the error correction procedure is repeated 
frequently enough. This can be demonstrated by the following explicit calculation of the state fidelity after error 
correction. 

In QECCs that correct single-qubit errors, the error operators are represented by of with I = 1,2, and 
a = 1,2,3 [11]. Let i? = I, denoting the identity operator, and i?3(z_i)+ Q = erf. We only consider orthogonal 
QECCs. Most of the discovered QECCs belong to this class [5-12]. For orthogonal QECCs, the encoded input state 
|* (0)) should satisfy the condition [10] 

(* (0)| Rl R n , |* (0)) = S nn , (n = 0, 1, • • • , 3L) . (17) 

During the error correction procedure, we first detect the error syndrome. If there is a R n error, i.e., the state becomes 
R n |* (0)), we apply the recovery operator R^ 1 = R n to the state and thus get the correct initial state |* (0)). In 
the case of spatially correlated decoherence, the error operators are represented by Q n (n = 0, 1, ■ ■ ■ , 3L), but we still 
adopt the above error correction procedure. If there is a Q n error, which occurs with probability p n , we detect the 
error syndrome and with probability |(\& (0)| R} n Q n |* (0))| 2 find that the error is R m . After this detection, the state 
is collapsed into R m \ ty (0)). We apply the recovery operator and thus get the initial state. The whole error correction 
procedure described above yields the following average state fidelity after error correction 

3L 3L 

F a (At) =^^p„|(f(0ptO B |f(0))| 2 (18) 

n—Qm—0 
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= 1-0 (At 2 ) . 

In deriving Eq. (18), we have used Eqs. (11) and (17), together with the identity \^k"\ 2 = 1 (from the unitarity 

l, a 

of the matrix U ). After the whole time T , the final average state fidelity F a (T ) is then approximated by 

F a (T ) ~ [F a (At)f ~l- (NAt 2 ) . (19) 

Since NAt 2 oc jj, it can be made arbitrarily small by a frequent repetition of the error correction procedure. This 
demonstrates that the QECCs devised to correct single-qubit errors are applicable for reducing spatially correlated 
decoherence. 

Before ending the paper, we should emphasize that we have omitted all the other terms except the first one in 
the Hamiltonian (1). If these omitted terms become important due to some special reason, the QECCs that correct 
for single-qubit errors may not work well any more. For example, if the second term in the Hamiltonian (1) is nor 
negligible, the quantum jump errors will include not only the terms that can be expressed as sums of single-qubit 
operators, but also the terms that are tensor products of them. To combat this kind of decoherence, the QECCs need 
at least having the ability to correct two-qubit errors. 
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